Physical model of continuous two-qubit parity measurement in a cavity-QED network 
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We propose and analyze a physical implementation of two-qubit parity measurements as required 
for continuous error correction, assuming a setup in which the individual qubits are strongly coupled 
to separate optical cavities. A single optical probe beam scatters sequentially from the two cavities 
and the continuous parity measurement is realized via fixed quadrature homodyne photo-detection. 
We present models based on quantum stochastic differential equations (QSDE's) for both an ideal 
continuous parity measurement and our proposed cavity quantum electrodynamics (cavity QED) 
implementation; a recent adiabatic elimination theorem for QSDE's is used to assert strong conver- 
gence of the latter to the former in an appropriate limit of physical parameters. Performance of the 
cavity QED scheme is studied via numerical simulation with experimentally realistic parameters. 

PACS numbers: 03.67.Pp,42.50.Lc,42.50.Pq,03.65.Yz 



It is now well established ^ that error correc- 
tion/avoidance protocols and fault-tolerant architectures 
are essential for any practical implementation of quan- 
tum information processing. While most theoretical re- 
search in these areas utilizes discrete time map-based 
models for quantum dynamics and decoherence, which 
are perhaps more familiar to computer scientists and in- 
formation theorists, there has been growing interest in 
transferring key ideas [2] to the domain of continuous 
time differential-equation-based models, which are more 
common in the context of ah initio physical modeling. 
In this paper we contribute to a line of research, initi- 
ated by Ahn and co-workers [EH, S 0] and broadened 
by other research groups 0, which focuses on con- 

tinuous quantum error correction via stabilizer coding 
and continuous syndrome measurement. This approach 
is attractive for design and analysis because it fundamen- 
tally connects the goal of quantum decoherence suppres- 
sion with formal optimization methods of classical con- 
trol theory 0]. It also has a significant potential imple- 
mentation advantage over standard discrete-time formu- 
lations in that continuous tracking of errors may be real- 
ized without the need for executing cumbersome readout 
circuits, but this of course relies on the assumption that 
continuous non-demolition syndrome measurement can 
be realized in an experimentally favorable way. In what 
follows we describe a straightforward implementation of 
continuous two-qubit parity measurement (sufficient for 
syndrome measurement of the quantum bit-flip code) in 
the context of cavity quantum electrodynamics (cavity 
QED), and analyze the performance of our scheme both 
for fixed realistic parameters (via numerical simulation) 
and in an ideal limit of parameter values (via adiabatic 
elimination). Our scheme utilizes a simple coherent-state 
optical probe in place of the usual ancillary qubits, and 
exploits Hamiltonian qubit-cavity couplings in place of 
clocked quantum logic gates for the syndrome readout. 
The strength of the syndrome measurement can never- 
theless be modulated easily (or even turned off entirely) 
by adjustment of the power of the optical probe beam. 



USA 




FIG. 1: Schematic depiction of two cavities driven sequen- 
tially by a resonant laser beam. A three level atom is trapped 
inside each cavity, and identical atom-cavity dynamics apply 
in each. After probing both cavities the laser light is directed 
to a homodyne receiver. 

The basic setup of our proposed implementation is 
shown in Fig. [TJ two optical cavities, each containing 
a single three-level 'atom' (potentially a gas-phase alkali 
atom, nitrogen vacancy center in diamond, etc.), are in- 
terrogated sequentially by a coherent optical probe with 
amplitude a (similar arrangements have previously been 
considered in the context of quantum information sci- 
ence [ioj). A qubit is encoded in the ground states | — ) 
and 1+) of the intracavity atom; an optical transition be- 
tween 1+) and the excited state |e) is coupled strongly to 
a quantized cavity mode with vacuum Rabi frequency g. 
For simplicity we assume atomic selection rules such that 
|e) decays only to |+), with excited state decay rate 7^, 
and single-sided cavities with photon decay rate . A 
homodyne receiver is used to measure the phase quadra- 
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ture amplitude of the laser beam after it probes the two 
cavities. The basic intuition behind our scheme is that 
the coherent probe acquires a phase shift of either tt or 
radians upon reflection from each cavity, depending on 
whether the intracavity atom is in the coupled or uncou- 
pled qubit state [11]. After reflecting from both cavities 
the probe carries an overall phase shift of tt radians if 
the parity of the two qubits is odd, and or 27r radians 
if the parity is even. As the latter two conditions are in- 
distinguishable, the homodyne measurement effectively 
implements a parity measurement. 

We model our system quantitatively using quantum 
stochastic differential equations (QSDE's) and quantum 
filtering (quantum trajectory) theory [l2[. The time evo- 
lution of an observable S of the atoms and fields is given 
by jt{S) = UlSUu where Ut = ^ dUs and 

dUt = 1^ + K:b2 + ot) dA\ — {nhi + t<ib2 + o^)^ dAt 

2 

i=i ^ 

+ (— (6i + ^2) - — (^1 + ^2)^) dt VUu (1) 

with f/o = This propagator acts on \)®^ T®'^ , where 
() = (8)^^(N) is the Hilbert space of a three- level atom 
and quantized cavity mode, and T is the bosonic Fock 
space of the probe and of the radiation modes for atomic 
spontaneous emission, hi and Gi act as the i^^ cavity 
mode annihilation operator and |+)(e| atomic lowering 
operator, respectively, and as the identity on the remain- 
ing spaces. Operators on the bosonic Fock spaces are 
represented as quantum noises At and B^^^ for the probe 
channel and i — {1,2} radiation modes. These annihila- 
tion processes are related to the more familiar Bose fields 
by At = J^asds, b[^^ = j^h^s^ds. 

Eq. ^ describes the evolution of the cavity QED sys- 
tem at its fundamental level, i.e., at the level of abstrac- 
tion often used in experimental considerations. Many of 
the more intrinsically scalable, solid state cavity QED 
implementations have additional dynamics not in ([1]), 
but which could easily be incorporated via additional 
dephasing terms. Although complete, ([T]) is unwieldy 
with operator coefficients that couple together an in- 
finite number of dimensions. In the strong coupling 
limit {y/g ^ n^a ^ 7) the effective dimensionality of 
the system is limited to certain slow degrees of free- 
dom and the fast dynamics can be adiabatically elim- 
inated from the description. This is accomplished in 
our formalism using the theorem of singular perturba- 
tions on QSDE's [l^, [ll|, from which it can be shown 



that the physical dynamics ([T]) is approximated (with 
Iii2 = Z1Z2, At = Jq alagds and Uq = Q ^ Ijrm) by 

dUt = {{U12 - Q)dAt ^ aUi2dAl - aQdAt 
lap 

-l^Qdt}Ut, (2) 

in the sense that V?/^ G Q J^^^, 

lim lim \\{Ut-Ut)^ =0. (3) 
f =const. 

Here Q projects onto Q, Q = Q^^, Qo = span{|ii), \d)}, 
\u) = 1+) |0), \d) = |-) where the cavity mode 

is here represented in a coherent amplitude basis, and Zi 
acts as \u){u\ — \d){d\ on the i^^ cavity system. Thus, the 
idealized dynamics represented by Eq. ([2j) live on a two- 
qubit Hilbert space coupled only to the probe. Note also 
that the qubit operators now appear as non-trivial three 
body operators simultaneously coupling both qubits with 
probe excitations (such as ZiZ2dAt). We have derived 
such effective, multi-body, non-local, qubit parity inter- 
actions as a very good approximation of the fundamental 
picture ([1]) in the strong coupling limit. 

Given a QSDE model, the conditional evolution equa- 
tions induced by continuous measurement can be de- 
rived straightforwardly [l2^. The resulting quantum fil- 
tering equation, also known as a Stochastic Schrodinger 
or Stochastic Master Equation, can be used to propagate 
a conditional quantum state that represents a sufficient 
statistic for the observer's best estimates of observables 
jt{X) on the basis of continuous measurement records 

If the conditional state remains pure for all times 
we may represent it by a vector G ()*^^. As complete 
purity of the conditional state can only be achieved if 
all output channels corresponding to Markovian environ- 
mental couplings are monitored with perfect efficiency, 
this is not a realistic assumption for laboratory imple- 
mentation. However the equations so derived are very 
useful for simulation and analysis; analogous equations 
for imperfect /incomplete observation with a mixed con- 
ditional state (density operator) are easily derived as re- 
quired for more practical purposes. 

If we assume homodyne detection of the quadrature 
amplitude observable f/^^ = jt{At + aJ) and photon 
counting of the radiation modes for atomic spontaneous 
emission Y^^"^ = jt{A^^^^) k = {1,2}, the 'physical' pure- 
state propagator (for the information state Vt) derived 
from ^ is 

dvt = l{f^h^i^b2^a)dY}^^ ^g^i^a^'^'^bi-a^'h]^^ 
[ i=i 
^2 

-{b\bi + b\b2)—dt - blbiK^dt - anibi + b2)Ut 

+ E((7.«-/)^r«-^^^./t)|.,. (4) 
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FIG. 2: Two independent simulations of atomic Bell- state 
preparation in a realistic system (see text). The top pair of 
graphs show a simulation in which the system happened to 
reduce to even parity; the bottom pair shows an odd-parity 
example. 



For simulation purposes Eq. (|4|) can be used to generate 
realistic homodyne measurement records dY^^^ by driv- 
ing the innovations and gauge processes with appropriate 
pseudo-random numbers, as is commonly done in quan- 
tum optics and atomic physics [15]. Similarly, if we define 
jt{X) = UlXUt using ([2]) we can construct the 'idealized' 
pure state propagator assuming homodyne measurement 
oiYt = ~jt{At^A\) as 



dvt = I aUi2dYt - ^-^Qdt 



vo e Q. (5) 



It is straightforward to show that the idealized filter ([5]) 
represents a finite-time unraveling of an ideal pro- 
jective measurement as discussed in [l6|. We can alter- 
natively think of ([5]) as a reduced filter for analyzing the 
homodyne photocurrent in our two-cavity setup, which 
exploits adiabatic elimination for a reduction in variable 
count. In particular we can track the parity of the two 
qubits (approximately, but with little computational ef- 
fort) by driving ([5]) with a given homodyne photocurrent 
in place of dYt . In an experimental scenario we would use 
the measured photocurrent; below we will also examine 
the results of driving ([5]) with simulated photocurrents 
generated by the physical Eq. (|4|). 

To demonstrate that our physical setup indeed re- 
alizes an approximate parity measurement, we use ([4|) 
to simulate entanglement generation from separable ini- 
tial states of the atom-cavity systems. We use param- 
eters that should be achievable in a Fabry-Perot cav- 
ity/cold Cs atom system with mm-scale dielectric mir- 
rors, {g,n^/2,-i'^/2,a} = {20,4.5,0.5,0.2}. We numeri- 
cally integrate (|4]) from the initial state '^o = 2~^((|+) + 
| — )) |0))'^^. As can be seen in Fig. [21 the system be- 
gins with equal projections on the 2~^/^(| + ±) + | — ^)) 
joint atomic states and the expected homodyne signal 
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FIG. 3: Summary of 1000 Bell-state projection simulations. 
The top graph compares the statistics of the ideal and physical 
parity variance. The mean of the Var(cr^^V^^'*) trajectories is 
shown as a solid, black line; grey shading indicates the one- 
standard deviation range of trajectories above and below the 
mean. The mean and standard deviation range of Var(ni2) 
are shown by the dotted line and blue shading. The bot- 
tom graph highlights the performance of the reduced filter. 
Referenced to the right axis, a representative Var(cr^^V^^'*), 
Varr/(ni2) trajectory pair are shown in solid black and dot- 
ted green, respectively. The blue-grey line and shading, refer- 
enced to the left axis, depicts the mean and standard devia- 
tion range of the fractional residual error over all trajectories. 



vanishes. The laser is then adiabatically switched on and 
the conditional state gradually projects into an atomic 
parity subspace with the corresponding (c/f/^^) ^ ±2a. 
After a fixed time j'^t = 80 the probe is adiabatically 
switched off. While the probe is on, the non-zero overlap 
with atomic Bell states not present in the initial state is 
caused partially by entanglement with the cavity mode 
states, but the residual expectation of such Bell states 
after the laser is switched off is wholly due to the ac- 
cumulation of weak measurements performed within the 
atomic parity subspaces by the field modes (see below). 

We now compare (j4]) and ([5]) by two different com- 
putational procedures. First, we compare indepen- 
dent simulations of the idealized projective measure- 
ment represented by Eq. ([5]) and the approximate pro- 
jections represented by Eq. (|4]). To this end, we con- 
struct trajectories of the variance of parity operators: 
Var((j^^V^^) from Vt in simulations of Eq. ^ (with 
crz = |e)(e| + |+)(+| - |-)(-| distinguishing cavity- 
coupled and -uncoupled states), and Var(ni2) from Vt 
in independent simulations of Eq. ([5|). Some summary 
statistics from these simulations are depicted in Fig. [3l 
The integrations of the idealized filter are initialized 
with Vo = 2-^{\u) + \d))^^. At t = we begin with 

Var(ni2) = 1, and both variances de- 
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crease in time as the systems randomly project into one 
parity subspace or the other. Indicative of their similar 
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statistics, the Var(cr^^V^^^) and Var(ni2) trajectory en- 
sembles largely overlap at all times. Moreover, it can be 
shown that the excited state population remains small 
at all times, (^)^, and that the atomic dynamics are 
principally constrained to the two ground states. 

Second, we assess how well (j5j) performs as a re- 
duced filter for analyzing the physical system. In this 
case, we construct trajectories of Varr/(ni2) by integrat- 
ing vt with photocurrents simulated from Eq. (|4]). A 
representative Var^j(ni2),Var(cr^^V^^) trajectory pair 
in Fig. [3] suggests the accuracy of the reduced par- 
ity estimate. Underlaying this are the statistics of 
the fractional residual error from 1000 such pairs: 
|Var(4^V^^VVar^/(ni2)|/Var(4^V^^^). Although the 
ranse of Var(4'Vz^) 

spans 5 orders of magnitude, the 
reduced filter performs well, tracking this physical parity 
estimate to within a factor of 2 in every shot. 

Now consider Var((j^^V^^) [ax = |-)(+| + |+)(-|), 
which serves as a measure of the distinguishability 
of the atomic states within each parity subspace: 
Var (cr^V^^) > indicates an imperfectly prepared 
atomic Bell state. The idealized description (|5|) makes 
no measurement within each parity subspace. Vt^ how- 
ever, is capable of (weakly) distinguishing states within 
parity subspaces through many different mechanisms. A 
lower bound Var(<j^'' cr^'') > 1 — exp(— 8a^//^^) is set by 
the non-separability of the atomic and cavity states in 
equilibrium. As the probe intensity is switched off, the 
systems factorize again. Spontaneous emission events, 
although rare (in each time step they occur with prob- 
ability ~ (^^)^(it), completely destroy any Bell state 
entanglement. But even if there are no spontaneous 
emission events the probe alone can make weak mea- 
surements within the parity subspaces in steady state. 
For exactly identical atom-cavities, the even parity sub- 
space has a weak constant decoherence mechanism: the 



expected homodyne increment dY^ of the l^)"^^ state 
is 2adt^ but is slightly less than that for the \u)^'^ state 
because of expected atomic scattering. If the cavities are 
estimated to be in the state cos(6>)|ii)^^ +sin(6>)|(i)^^, the 

^-information content of the next dY^^^ can be quantified 
by its Fisher information I{e) = {2as\n(2e){'^ffdt. 
The larger T{0) is, the more likely our estimate of 9 will 

change because of dY^^\ Note that T{6) is maximized 
for the Bell states in this parity subspace, making them 
especially fragile even in steady state. Of course identi- 
cal cavity systems are an idealization, but while asymme- 
tries between the cavities can decrease the fragility of the 
even-parity subspace they increase that of the odd-parity 
subspace. In general, probe fluctuations cause decoher- 
ence as all I ± ±) states respond differently to perturba- 
tions. Symmetric-cavity considerations of Var(cr^V^^) 
are thus useful as best- and worst-case scenarios of steady 
state, probe-induced decoherence. 

In summary, the close approximation of ideal parity 
measurement by a setup based on coherent probe fields 
and homodyne detection establishes the two cavity sys- 
tem as a continuous-time alternative to the common 
quantum-circuit building block of sequential controlled- 
not gates with an ancillary qubit. In this paper we have 
demonstrated the use of an adiabatic elimination the- 
orem for QSDE's and numerical simulation of derived 
quantum filtering equations to assess the performance of 
our proposed implementation of continuous parity mea- 
surement in both ideal and realistic parameter regimes. 
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